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MEASURING THE MIXING EFFICIENCY IN A 
SIMPLE MODEL OF STIRRING:SOME ANALYTICAL 
RESULTS AND A QUANTITATIVE STUDY VIA 
FREQUENCY MAP ANALYSIS 

TIMOTEO CARLETTI AND ALESSANDRO MARGHERI 


Abstract. We prove the existence of invariant curves for a T- 
periodic Hamiltonian system which models a fluid stirring in a 
cylindrical tank, when T is small and the assigned stirring protocol 
is piecewise constant. Furthermore, using the Numerical Analysis 
of the Fundamental Frequency of Laskar, we investigate numeri¬ 
cally the break down of invariant curves as T increases and we give 
a quantitative estimate of the efficiency of the mixing. 


1. Introduction 


In pO it is studied a simplified model of the stirring of an ideal fluid 
in a cylindrical tank by an agitator. A Lagrangian representation is 
considered for the motion, which is assumed to be completely two di¬ 
mensional. The tank thus degenerates to its boundary circle, which has 
radius R. The agitator is modelled as a point vortex of strength T and 
its position inside the boundary circle as a function of time, denoted 
by z(t), is a prescribed T-periodic function called stirring protocol. In¬ 
troducing the complex coordinate ( = x + iy, the motion of the fluid 
is governed by the following nonautonomous T-periodic Hamiltonian 
system 


( 1 ) 


t r \ z (t )\ 2 — R 2 


whose corresponding Hamilton function is given by 


( 2 ) 





C - z(t) 
z{t)C - R 2 


System (JTJ) is defined on the set {(£, () e RxC : ( ^ z(t),( 7^ = } 
and it is considered for ( belonging to the invariant disk Dr = {( 6 
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C : |C| < -R}. The following piecewise constant stirring protocol 


( 3 ) 


z{t) 



t £ nT, nT + T/2^j 

t£ nT + T/2, (n + 1)T^ n £ TL , 


with fixed 0 < b < R, has been investigated in JTJ. In this case the 
motion can be integrated over finite time, whereas the long time be¬ 
havior of the model is studied by numerical experiments when 6=1/2 
for different values of T. In JTJ, the objective of the author is to study 
which features control the onset of chaos and hence the efficiency of 
the mixing. For this aim, the regimes of regular and chaotic behav¬ 
ior are identified and a qualitative description (see JT] page 736 second 
paragraph) of the mixing efficiency of the model has been provided for 
several values of the parameters ( b,T ). It turns out that the regions 
of chaotic behavior consume a larger and larger portion of the phase 
space as T increases, disrupting the regular pattern observed when T 
is very small (see Figure QJ. Quoting JTJ, '...for small T the model [...] 
will look more and more like the two fixed agitator system; [...] Thus 
one would expect ’convergence’ as T 0. ’ 

In Section 2 of this work we will show that averaging theory pro¬ 
vides the appropriate framework to investigate the above statement. 
Loosely speaking, in such a setting we will be able to show that when 
T —> 0 the flow of system 0 with the stirring protocol © converges 
in the C A topology to the Hamiltonian flow corresponding to two fixed 
agitators. As a consequence, by using Moser’s Small Twist Theorem 
(|E|,|I1.|7|,0) we will prove that the regularity of the pattern ob¬ 
served in the numerical experiments for small T (see Figure 2 a), b) 
and c) page 733 of |T] and Figure d for T = 0.05) is due, at least in a 
suitable annular region inside Dr, to the presence of invariant curves. 
In fact, such curves are, as one readily realizes, an obstruction to the 
mixing of the fluid. 

The result which we present with this respect is the following: 



FIGURE 1. Different behaviors of the Poincare map, 
with b = 1/2, R = 1 and T = 0.05 (left), T = 0.5 (mid¬ 
dle) and T = 1.5 (right). Each plot is obtained following 
10000 time units the same 18 initial data. 

Theorem 1. Let z(t) be the stirring protocol ©) mid let T > 0 be its 
period. Then, there exists an annular region A inside Dr such that the 
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T-Poincare map associated to system m has invariant curves in A for 
any sufficiently s?nall T. 

As T increases the invariant curves break down. In Section 3, we will 
study numerically this phenomenon by means of the Numerical Anal¬ 
ysis of the Fundamental Frequency, NAFF of Laskar HU- This allows 
us to extend our previous analytical result to cover larger T ranges, by 
showing numerical evidence that invariant curves persist close enough 
to the boundary of the invariant disk Dr for large values of T. Further¬ 
more, we are able to give a quantitative description of the efficiency of 
the mixing in function of the parameters (b, T ), by measuring, through 
NAFF, the portion of phase space filled by invariant curves. 

Acknowledgements. Support from GRICES/CNR project is acknowl¬ 
edged. The second author would also thank support from FCT 

2. Existence of invariant curves 

In order to establish the existence of invariant curves for system 0 
with the piecewise constant protocol © and small T, we will show that 
its T-Poincare operator (i.e. its period map) satisfies the assumptions 
of the so called Moser’s Small Twist Theorem. 

Before outlining the strategy of our proof, it is worth to recall the 
statement of this important result. 

Fixed 1 1 > Iq > 0, let us consider the annulus defined by A = 
{(0,1) : 0 G S 1 , 1 0 < I < Ii}. Given a mapping M : A —> M 2 , we can 
End a lift to the universal cover A — {(0,1) : 0 G M, J 0 < / < ii} of 
A, which we still denote by M. 

Moser’s Small Twist Theorem Let a be a G 4 ([/ 0 ,/ 1 ]) function sat¬ 
isfying 

(4) «'(/)< 0 V/ g [J 0 , h]. 

Then, there exists e > 0, depending on I\ — Iq and a, such that the map 
M : A ¥ M 2 has invariant curves if it satisfies the conditions below 

a) M has the intersection property, that is, for any Jordan curve 
7 homotopic to the circle I — I 0 in A, M(^y) fl 7 ^ 0/ 

b) the lift of M can be expressed in the form 

(5) M(0,1) = (0 + Ta{I) + TMO , I,T),I + Tfc(0,1, T )), 

for some T 6 (0,1) and £ CS(A) with ||</>i||c 4 (A) + 

1102||c 4 (A) < e - 

This version of the theorem is presented in na and may be proved 
using the techniques developed in urn- 

In what follows, a map of the form (6 + Ta(I),I) with a'(I) 0, 

/ G [/ 0 ,/i], will be referred to as small twist map. 

The application of the previous theorem will be obtained perform¬ 
ing the following steps. First, by using a classical construction from 
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averaging theory, in Subsection 12.11 we will rewrite system m as a 1- 
periodic perturbation of the Hamiltonian averaged system, being T the 
small parameter. The averaged system corresponds to the two fixed ag¬ 
itator model. Then, exploiting the geometry of the phase space of the 
averaged system, we will construct explicitly the corresponding action 
variable on a suitable subset of Dr. As a consequence, we will be able 
to show that on an appropriate annulus inside Dr the 1-Poincare map 
of the averaged system is a small twist map. 

Next, the C 4 estimates needed to apply the Small Twist Theorem 
will be provided by a general result about the differentiability of a 
flow with respect to parameters, which we recall in Proposition [21 of 
Subsection 12.21 

Finally, in Subsection 12 .Ml we collect all these fact and complete the 
proof of Theorem 1. 


2.1. The averaged system. In what follows, we identify C with M 2 
and ( — x + iy with (x,y). We rewrite the Hamilton equations of 
system © in a compact real form as 


( 6 ) C = JV ( H(t,C), 

where J is the standard 2x2 symplectic matrix and J^). 

The Hamilton function H (t, Q is piecewise autonomous and may be 
expressed in the form: 

Hit, o = mn+i o + (i - m)H-i o , 


where 

(7) 


H±{C) 


ill CTfr 

27T 11 b( =f R 2 


and 0(f) is the T-periodic extension of the restriction to [0, T] of 
the characteristic function of [0,T/2], Hence, defining V := M 2 \ 
{±6, ±R 2 /b}, H(t,() is smooth on the set (M \ {kT/2,k G Z}) x V. 
Moreover, H and its derivatives with respect to which exist in V for 
any t G M, have jump discontinuities at f = kT/2, k G Z. 

Next step is to consider T as a small parameter and to determine a 
suitable comparison limit system when T —> 0. For this aim, as in [T|, 
we first rescale time in © by f = Tr, so normalizing to 1 the period of 
the stirring protocol and of the vector field. In the new time r, setting 
Hi(C,t) := H{(,tT) and y(r) := C( r 70, system © takes the form 


(8) y(r) — TJVyHiir, y). 

To this system the theory of averaging applies (see [B] [Section 4.4]) as 
follows. Denote by 

(9) Hxiy) ■= H x (r, y) dr = ^ In 


y 2 — b |2 
y 2 - R 4 /b 2 
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and let V C C be an open set such that its closure is contained in 
D. Then for small enough T 0 > 0 we can find a symplectic 1 , close to 
identity, change of coordinates of the form 

(10) y = ij + Tw(t,ij), (t, 77, T) G M x V x (0, T 0 ], 

with w 1 -periodic in r and such that 


(11) w(0,77) = 0, 

which transforms system (JHJ into system 

(12) f] = TJV r] Hi(r]) + T 2 h{r, rj, T). 

By construction, the function h(r,r),T ) in (TH21) is 1-periodic in r, has 
jump discontinuities at r = fc/2, k G Z and it is smooth on (M \ 
{/c/2, k G Z}) x V x [0, T 0 ]. 

System (JBJ) is a perturbation of the following (integrablc) averaged 
system 

( 13 ) r } = TJV TI H 1 (v)- 


Henceforth, for simplicity, we will write H instead of H\. 

We note that the Hamilton function H corresponds to a two point 
vortices system, one vortex being located at (+ 6 , 0 ) and the second at 
(—b, 0). Moreover, exploiting the geometry of the phase space of system 
(USD, we can construct explicitly the action variable for this system 
outside the homoclinic loops surrounding the vortices (see Figure El on 
the left for T = 0.05). This is done as follows. 

By introducing symplectic polar coordinates ( = \[2re 2 ^ and setting 
for notational convenience E = e 47rH ^ r we can express the level lines of 
H outside the homoclinic loops 


(14) r(-0) = 


b 2 R A 
R A + 6 4 


cos 2-0, 0 G 


■~r 7r , t 71- 
4 4 


U 


—71, -7T 

4 4 


in the form r = r(0,T) with 0 G M and E G (b 2 /R 4 ,l/R 2 ]. We 
observe that E = 1 /R 2 corresponds the boundary of Dr, whereas 
b 2 /R 4 corresponds to the homoclinic loops. 

Taking into account that r(0, E) = r(—0, E) = r (tt — 0, E) (namely 
the averaged system is invariant with respect to ( —> ( and ( —> —0), 
the action variable / is given by 
(15) 


1(E) = 


J H(V2re i *)=^lnE 


r dip = 


<E) 

2vr 


fir/2 


a Jb(E) cos 2 (20) + c(E) dip , 


Tntroducing canonical variables y = ( p , g) and ?? = (P, Q) the requested trans¬ 
formation can be obtained trough the following generating function S(q,P,r) = 

Pq + T J 0 T Hi(P, q) — H\{P, q , s) d.s. Thus p = d q S can be inverted if T is small 
enough and trivially S(g, P, r + 1) = P(g, P, r) for all (q, P, r) which implies (till) . 
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where 

a(E) := 1 - W ' MC := 46 4 (1 - E 2 /? 1 ) 2 
c(£) := 4(£V - 1)(6 4 - E 2 R 8 ) . 


In next proposition, we use the action variable constructed above 
to show that the 1-Poincare operator of the averaged system, which 
henceforth we denote by Mf, is a small twist map on an appropriate 
annular region inside Dr. 


Proposition 1. The set 

„ . ^ , b 2 1 

f := Be 


: T(E) ± 0 and El" (E ) + T(E) ± 0 \ , 


A 4 ’ R 2 

is a non empty interval. Fix E 0 e (inf S, -^) and consider the annulus 


os) ' 4 :={c efl * : s ln£; “-" (c) -£ ln F}' 

Then, Mj is a small twist map in A. 

Proof. An easy computation shows that 1/R 2 G £. The rest of the 
statement about £ follows from the continuity of the maps E t —> I'(E), 
E i—> EI"(E) + I'(E). Let us consider action-angle coordinates (0,1) 
for the averaged Hamiltonian system (I13[) . with I given by (1T51) . In 
the annular domain A, defined by m, fi ie ma p i = -f(-fB) ; = i{e(h)) 
has an inverse H = H(I ) on the interval [J 0 , R 2 / 2], where Jo := I(E 0 ). 
Moreover, on such interval we have H"(I) ^ 0. Since, H”(l/R 2 ) < 0, 
it follows that 

(17) H"(I)< 0, Ie[I 0 ,R 2 /2}. 

The 1-Poincare map associated to system m admits a lift to the set 
A := Rx [I 0 , R 2 /2] given by Mf (6,1) = (0 + TH'(I), I) and therefore, 
taking also into account CED, it follows that it is a small twist map. □ 

In order to apply Moser’s Small Twist Theorem, we need some fur¬ 
ther estimates, which will be provided by the consequences of the 
Peano’s Theorem stated below. 

2.2. Differentiability with respect to parameters. Let T > 0 be 

the period of the stirring protocol © and let Mf be the 1-Poincare op¬ 
erator of system In this subsection we will recall a general result, 
namely Proposition El below, which is a consequence of the differentia¬ 
bility of the flow of a system of O.D.E.’s with respect to parameters. 
This proposition is a restatement of Proposition 6.4 in m under the 
standard assumptions considered in dealing with vector fields which 
are discontinuous in the independent variable. This result will provide 
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the estimates needed to show that, for sufficiently small T, the maps 
Mf and Mf are C 4 -close on the annulus A defined by m- 

This will be done in the next subsection, and the proof of our main 
result will follow easily. 

Let V be an open subset of R n , let T 0 be a hxed positive number 
and let m be a nonnegative integer. Consider the following differential 
equation depending on one parameter 

AT 

(18) - = F{t,x,T), 

where F : [0,1] x V x [0, T 0 ] —> M n , (t, x, T) > F(t, x, T ) is a function 
which satisfies 


• the map (x,T) i—> F(t,x,T ) £ C m+1 (V x [0, T 0 ]) for almost all 
*e[ 0 ,l]; 

• for 0<j + k<m + l, the maps t ^ af ^jF(t,x,T) are t 
measurable for any (x,T) £ V x [0,T 0 ]; 

• for each compact set K c V there exists C > 0 such that 


Qk+j 

dx k dTi 


F(t, x, T) 


<C, 


for all 0<k + j<m + l and (t, x, T ) £ [0,1] x K x [0, T 0 ]. 

The solution of m satisfying x(0) = x 0 , will be denoted by x(t] x 0 , T). 

By the general theory of ordinary differential equations, x is of class 
C 0,m+l,m+l 

in its three arguments whenever it is defined. The following 
result is a consequence of this fact. 


Proposition 2. Let A be a compact subset of V such that for every 
xq £ A and T £ [0, To] the solution x(t',xo,T) is well defined in [0,1]. 
Then, for each (t, xq, T ) £ [0,1] x A x [0, T 0 ] the expansion below holds 

dx 

x(t] X 0 , T ) = x(t] Xq, 0) + — (t; Xq, 0)T + R(t ; Xq, T)T , 
where the remainder R satisfies 

•>7 n )||c m (A) —>■ 0, T —■> 0 , 

uniformly in t £ [0,1]. 

2.3. Proof of Theorem ^ We are now able to prove our main the¬ 
orem. We will show that for sufficiently small T the T-Poincare map 
associated to system 0 with the stirring protocol © has invariant 
curves in the set A defined by GU). 

For a hxed R' satisfying R < R' < R 2 /b, we define 

, 7,2 734 , 

(19) V := |C = V2re^ £ C : R4 + &4 cos 2 if < r < (i?') 2 / 2 } • 

This set, which contains A, corresponds to the open disk of radius 
R! minus the compact region bounded by the homoclinic loops. Hence, 
it has positive distance from the singularities of H. Now we choose 
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a sufficiently small T 0 > 0 for which the transformation m is well 
defined on [0,1] x V x [0, T 0 ] and, moreover, all the solutions of system 
© with initial conditions in A are well defined in [0,1]. Let T < T 0 
be the period of z(t). By Proposition [21 and by introducing the action 
angle coordinates considered in Proposition El we can find a lift of Mj 
to A := R x [Jo, R 2 /2 ] of the form 

Ml (0, /) = [Q + H\I)T ,, /) + R(6, /, T)T , 

where ||/?(-,T)||c4(^) —> 0 for T —> 0. Hence, Ml satisfies hypothesis 
b) of Moser’s Small Twist Theorem for T small enough. Moreover, Mj 7 
has the intersection property in A, being an area-preserving map in A 
for which the boundary of A is invariant. 

Then, by Moser’s Small Twist Theorem, Ml has many invariant 
curves in A for T small enough, and so does the T-Poincare operator 
of system m, which obviously coincide with Mj. 

3. Numerical Results 

In this section we will study system m with the stirring protocol © 
from a numerical point of view. More precisely, we will be interested 
in two types of numerical experiments. In the first one, we will explore 
the continuation properties with respect to T of the invariant curves 
of the T-Poincare map of system © obtained for small T in Theorem 
□ In S3 we will give numerical evidence of their persistence close to 
the boundary for large values of T. 

In the second experiment, we will be able to give a quantitative de¬ 
scription of the mixing efficiency of the piecewise constant stirring pro¬ 
tocol. This will be done in H3.3I where we will investigate the parame¬ 
ters (6, T) giving rise to an efficient stirring. 

Both these numerical experiments will be carried out by using the 
Numerical Analysis of the Fundamental Frequency (NAFF) of Laskar[TU] 

This is a numerical method which allows to obtain a global view of 
the behavior of a dynamical system by studying the properties of the 
frequency map, numerically defined from the action-like variables to 
the frequency space using adapted Fourier techniques. This method has 
been used to investigate a wide class of dynamical systems, such as the 
solar system d, the galactic dynamics na, particle accelerators [T2] 
and the standard map 

For sake of completeness, let us first present briefly the main outlines 
of the method, referring to ra for further details. 

3.1. Frequency Map Analysis. Let us consider an n-degrees of free¬ 
dom quasi-integrable Hamiltonian system, expressed in action-angle 
variables by 

( 20 ) 


nn.n-,) = //, l i/i i/,«i, 
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where H is a real analytic function of (/, 9) G B x T n , B is an open 
domain in M n , T n is the n-dimensional torus and e is a small real 
parameter. For e = 0 the system is integrable: the motion takes place 
on invariant tori Ij = Ij{0 ) described at constant velocity Uj (/) = 

, for j = 1 ,,n. Assuming a non-degeneration condition on 


8H o 

a/, 


1 ( 0 ) 


H 0 , the frequency map F : B —> M n 

F :I^F(I) = v, 

is a diffeomorphism onto its image In this case KAM theory mm 
ensures that for sufficiently small values of |e|, there exists a Cantor 
set C of frequency vectors satisfying a Diophantine condition, for 
which the quasi-integrable system with Hamilton function stiii 
possess smooth invariant tori. These tori are e-close to those of the 
unperturbed system, and support the linear flow t > 9j ( t ) = Ujt + 
9j (0) ( mod 2n) for j = 1,... ,n. Moreover, according to Poschel 
there exists a diffeomorphism T : T n x —> T n x B 

*:(<£, i/)->(0,7) , 

which is analytic with respect to (j) in T n and C°° w.r.t v in fl e , and 
which transforms the Hamiltonian system generated by (EDI) into 

M = o 

(i) = Vj . 

For frequency vectors v G O e , the invariant torus can be represented in 
the complex variables (zj = IjF° :i )j = \ in by a quasi-periodic function 


( 21 ) 


z j (t) — z j (0) e J + a j,m iy) e 




Taking a section 9 = 9 0 of the phase space, for some 9 0 G T n , we obtain 
the frequency map Fq 0 : B Q 


( 22 ) 


Fe, ■■ 1 1 - n 2 (f- 1 (So, 0) , 


where 7 t 2 (</>, v) = v is the projection onto fh For sufficiently small |e| 
the non-degeneration condition ensures that Fg 0 is a smooth diffeomor- 
phism. 

If we have a numerical (complex) signal over a finite time span 
[—A", K} and we want to recover a quasi-periodic structure, we can con¬ 
struct an A^-terms quasi-periodic approximation / (t) = J2k =1 . 

Frequencies and amplitudes are determined with an iterative 
scheme (possibly) involving some weight function (Hanning filter). 

Assuming some ’’good” arithmetic properties of the frequencies and 
using the p-th Hanning filter, then one can prove jTUj that NAFF 
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converges towards the first 2 ’’true” frequency of the given signal, zq, 
with the following asymptotic expression for K —> +oo 

(») 

which will usually be several order of magnitude better than the O (1 /K) 
order obtained with FFT. 

Using the NAFF algorithm, it is possible to construct numerically a 
frequency map (see Figure EJ) in the following way: 

a) fix all angles to some value 6*o; 

b) for all initial values Jo of the action variables, integrate numeri¬ 
cally the trajectories with initial condition (Jo, #o) over the time 
span K: 

c) look for a quasi-periodic approximation of the trajectory, with 
the previous algorithm identifying the fundamental frequency 
zq^ of this quasi-periodic approximation. 

Because of (1231) . we can use NAFF algorithm to test the goodness 
of the reconstructed quasi-periodic approximation of the given signal, 
showing numerical evidence that motion doesn’t take place on invariant 
tori because of some ’’diffusion in frequency space”. More precisely, 
this can be done by replacing step c) of the previous algorithm with 
the following 

c’) divide the time span [—A', K\ into smaller parts (possibly over¬ 
lapping): [ti,ti + Ad], for some tj. Then NAFF reconstructs for 
each l a frequency iA Kl \l). If the frequencies (v^ Kl \l))i coin¬ 
cide, up to some prescribed numerical precision, when / varies, 
then we can conclude that we obtained a good quasi-periodic 
approximation. Otherwise, we can measure the diffusion in the 
frequency space. 

3.2. Application I: existence of invariant curves. Let us con¬ 
sider the problem of the existence of invariant curves homotopic to the 
boundary of Dr for the T-Poincare map of system m with stirring 
protocol ©• One readily realizes that every such curve will intersect 
the vertical segment £ := G C : x = 0, y G [0, A]}. 

In order to reconstruct a frequency map, we fix a piecewise constant 
stirring protocol by choosing the parameters (6, T). Then, for N initial 
data on the segment £ we construct numerically the T-Poincare map 
and we iterate it Ni ter times. 

Remark 1 (Numerical computation of the Poincare map). As already 
remarked in [OP, one can construct the Poincare map for the piecewise 
constant protocol avoiding the integration of the Hamiltonian system 

2 We assume to enumerate frequencies according to decreasing amplitudes: 
\a^\ > l a i+il- That is, some smoothness of the signal is assumed. 
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m, which results CPU-time consuming and introduces additional er¬ 
rors. The idea is that in each half period the Hamiltonian system is 
autonomous, hence integrable and, moreover, the integration can be ex¬ 
plicitly done: for t G [0, T/2), the motion takes place on an arc of 

circle with radius p = (y- ~ bj, centered at ( c = yyfyj , with 


C(0 )-b 


Thus, setting £(0) = ( c + pe l ^°\ after half period the 


^ ^ K(o)-R 2 

point will be at ((T/2) = ( c + pe^UU) where: 


2 A r 1 — A 2 T 

( 24 ) 

The motion in the second half period is similar. 

Equation \2J/\) can be solved using a sixth order Newton’s-type method. 
However, when b is close to R and/or T is large, Newton’s algorithm 
doesn’t converge to the good solution unless we provide a very accurate 
approximation of <f>(T/2) as initial value for the algorithm. In order to 
overcome this difficulty, we used the following trick. Equation is 
nothing but Kepler’s equation and in Id an explicit solution is given 
in terms of Bessel’s functions (which results, in computations, less ef¬ 
ficient than our Newton’s-type method). Neverthless for large values of 
the parameters we obtain a good approximation to start with Newton’s 
algorithm by summing few terms of the previous explicit solution. 


Once we have an orbit, (Mf((( 0))) 0<fc<jV , we use the NAFF al¬ 
gorithm to reconstruct a quasi-periodic approximation of it. The fre¬ 
quency map is not smooth at all because of the presence of resonances: 
there is no invariant curve with rational rotation number. Thus, for 
a fixed resolution, the existence of invariant curves is assumed if the 
frequency map looks ’’regular” (see Figure [2] for T = 0.05), whereas to 
’’irregular” graphs we associate non existence of invariant curves (see 
FigureOlfor T G {0.5,1.0,1.5}). 

In Figure [21 we present some numerical results for the piecewise con¬ 
stant protocol © with 5 = 1/2 and several values of T. On the left, for 
very small T, say 0.05, the curve looks regular, hence a large part of 
the disk is occupied by invariant curves. We can observe that the map 
is not twist in the whole disk: there is a point where the derivative of 
the frequency map is zero. For slightly larger T, say T G (0.125,0.2), 
the curve is still regular, but an irregular pattern is showed close to the 
origin. Roughly speaking, no invariant curves are present for T = 0.2 in 
the disk |£| <0.22. Also, elliptic points due to resonances are showed. 

In Figure 01 when T increases, the irregular patterns consume larger 
and larger portions of the invariant disk, but still a regular frequency 
map is present close to the boundary of D R . 
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Figure 2. The Frequency Map numerically recon¬ 
structed by NAFF. Parameters are 6 = 1/2, R = 1 
and T G {0.05,0.125,0.2}. Initial data are 1000 points 
(/•(O') = ir equally spaced on [0, i] and orbits are corn- 
put 



FIGURE 3. The Frequency Map numerically recon¬ 
structed by NAFF. Parameters are 6=1/2, R = 1 
and T G {0.2, 0.5,1.0,1.5}. Initial data are 1000 points 
CifO) = irj, equally spaced on [0, i] and orbits are com¬ 
puted for N iter = 50000. 

This behavior of the system is displayed in Figure 0 where we show 
a numerical result for 6=1/2 and T = 3.0, a period which is 60 times 
larger than the smallest value in Figure El Note the change in the scale 
on the axes w.r.t. to Figure El and 01 (the scale on the vertical axis has 
been changed to show the large oscillations of the frequency map). 

Other results for larger values of T exhibit a similar behavior. This 
supports the conjecture that invariant curves persist, close to the bound¬ 
ary, for arbitrarily large values of T. 








13 



Figure 4. The Frequency Map numerically recon¬ 
structed by NAFF. Parameters are b — 1/2, R — 1 and 
T = 3.0. Initial data are 1000 points C/(0) = irj, equally 
spaced on [0.98«, i] and orbits are computed for N iter = 
500000. 


3.3. Application II: mixing efficiency. We have already observed 
that invariant curves are an obstruction to global mixing. We note now 
that close to robust invariant curves (i.e. associate to good frequencies) 
there is a neighborhood filled by invariant curves El El, which prevents 
also from local mixing. Using the precision of NAFF and observation 
c’) in page m we can test the efficiency of the stirring protocol by 
evaluating the goodness of the reconstructed signal. 

More precisely, we choose a couple of parameters ( b,T ), then we di¬ 
vide the invariant disk Dr into a fine grid of N tot points. For each point 
we determine numerically (see Remark H) the orbit for a time interval 
[0,/l], Then, fixing some K± < K and some positive (f;); we divide 
the orbit into pieces (C ti+K{\ suc ^ that intervals [/;,/; +Ad] over¬ 
lap. On each piece, the use of NAFF gives us a fundamental frequency 
z/( A d(/), where we emphasize the dependence of the frequency on the 
/-piece of orbit. 

For a given orbit with initial datum £ 0 = (x 0 ,y 0 ), let us define 


(25) 


e(xo,yo) 


-log 


^ max; v^\l) — min; 

max; zWW(Z) + min; v( Kl \l) 


This is a good indicator of the robustness of the orbit. In fact, if 
£ (x 0 ,y 0 ) is ’’large” then max; ^ Ai )(l) and min; v^ Kl \l) are ’’close to¬ 
gether”. Thus, the reconstructed frequencies on different pieces of 
orbit do not vary too much, and we can assume that we are on a 
quasi-periodic orbit. On the other hand, if e(x 0 ,y 0 ) is ’’small” then 
max; and min; u^ Kl \V) are ’’far from each other”. In this case, 

since the reconstructed frequency is not constant, we conclude that we 
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are not on a quasi periodic orbit. Intermediate values of e(xo,yo) give 
rise to intermediate degrees of robustness. 

We now fix a threshold value £ t hr > 0 and we measure the portion 
of initial data in Dr to which one can associate a e^-robust orbit 

#{(>0)?/o) : e(x Q , Vo) > Zthr} 

A preliminary analysis, see Figure 0 of some robust orbits allows us to 
choose an appropriate threshold value. 





/ / 


FIGURE 5. Distance between the frequency v[ Kl \l), de¬ 
termined by NAFF, and the mean value of the u[ Kl \iy s, 
in function of the part of orbit [ti,ti + K\\. On the left 
very regular orbits, on the right less regular ones. Both 
plots are made with N iter = 100000, K\ = 50000 and 
ti = 50(Z - 1). 


In this way we are able to turn the qualitative analysis of Figure 7 page 
738 of [I] into a quantitative one. 

With this respect, we present two results, the first with e thr = 12, cor¬ 
responding to max; — 1 


min; V 


and 


< T i Hi) 

< 1.510 


-7 


< 3 10 6 , and the second with e thr = 15 
Each result has been obtained divicl- 


max; i/( T i) 
min; i4 T i)(/) 

ing the invariant disk in approximately 30000 points, equally spaced 
in both x and y by 0.01. Then, several values of b and T have been 
considered. 

Finally, by using the previous ideas, we precise as follows the classi¬ 
fication scheme given for a regime in [I]: 


[I] integrable if 0.6 < m £thr < 1 (very poor mixing property) ; 
[T] transitional if 0.3 < m £th < 0.6; 

[C] chaotic if 0 < m ethr < 0.3 (efficient mixing), 

and we summarize our results in the following Figure |b] 
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Figure 6 . A part of parameter plane and the corre¬ 
sponding m £thr function; on the left mi 2 , on the right 
mi 5 . Squares correspond to Integrable case, stars to 
Transitional whereas circles to chaotic. 

4. Conclusions 

In this paper we proved that the simple stirring model given by © 
has invariant curves for all T small enough and every b G (0,i?); the 
use of NAFF gives us numerical evidence that such invariant curves 
persists even for T large, closer and closer to the boundary of the disk 
as T increases. Hence, for small T the mixing efficiency in not good at 
all, whereas for larger values of T, say T > 1, it becomes reasonably 
good once b > R/2. Thus, even if very simple, the model can give rise 
to efficient mixing. 

We note that, due to the general results on which it relies, our proof 
of the existence of invariant curves for small T may work for more 
general T-periodic stirring protocols z(t). However, the starting point 
should be to obtain a ’simple’ averaged system. 

From the numerical point of view, the more general systems referred 
above can be surely studied by means of NAFF method, once we ob¬ 
tained good integrators for these systems. Still considering our piece- 
wise constant model, we think that using NAFF one could give precise 
estimates of the size of the annular domain close to the boundary of 
Dji containing invariant curves and of the rate at which it shrinks to 
zero when T increases to infinity. 
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